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I.  INTRODUCTION 


A  maximum  flow  network  is  defined  hv  a  set  of  arcs  and  a  set  ot 

points  called  nodes.  Each  arc  joins  two  nodes  and  has  associated  with 

it  a  positive  capacity  which  represents  the  maximum  amount  of  flow 

that  may  pass  over  it.  One  of  the  nodes  is  designated  as  the  source 

and  another  as  the  sink.  From  these  nodes,  arcs,  and  capacities  the 

maximum  amount  of  flow  that  may  pass  from  source  to  sink  may  be 
•k-k 

calculated. 

This  investigation  is  concerned  with  a  sensitivity  analysis  on  a 
class  of  such  networks  knowm  as  planar  networks.  Specifically,  each 
arc  of  the  network  is  subject  to  anywhere  from  one  to  n  breakdowns 
which  result  in  a  reduction  in  its  capacity.  The  amount  of  this 
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reduction  in  capacity  is  a  random  variable  with  known  mean  and  variance. 
It  is  desired  to  find  the  smallest  possible  value  of  F  resulting  from, 
at  mostyn  breakdowns  where  F  satisfies: 

P{ma::  flow  it  F}  £  B 

where  P  stands  for  probability  and  B  is  a  small  number  between  zero  and 
one . 


II.  MAX- FLOW  MIN-CUT  THEOREM 


The  central  idea  in  the  solving  of  maximum  flow  network  problems 
is  summarized  in  the  max-flow  min-cut  theorem.  First,  it  is  necessary 
to  define  a  cut  set  and  its  value. 

De f in i t ion :  Consider  a  network  consisting  of  nodes  that  include 
a  source  and  a  sink,  and  capacitated  arcs  that  join  two  nodes.  Let 
A  and  B  be  a  partition  of  the  nodes  such  that  the  source  is  in  A  ami 
the  sink  is  in  B.  Then  the  set  of  arcs  that  join  a  node  in  A  to  a  node 
in  B  is  called  a  cut  set  and  is  denoted  Ca.b]  .  Furthermore,  the  value 
of  this  cut  set,  V[A,B]  ,  is  equal  to  the  sum  of  the  capacities  i-f  its 
a  res  . 

According  to  the  max- flow  min-cut  theorem,  the  maximum  flow  equals 

rsi 


the  minimum  value  of  all  cut  sets. 


III.  THE  TOPOLOGICAL  DUAL 


The  topological  dual  of  a  network,  when  defined,  is  another 
network  in  which  the  arcs,  instead  of  having  capacities,  have  lengths. 
Furthermore,  there  is  a  one-to-one  cor  resDondence  between  the  proper 
cuts  of  the  original  network  and  the  routes  through  the  dual,  and  the 
problem  of  finding  the  minimum  cut  may  be  reduced  to  one  of  finding  a 
shortest  route. 

Let  the  original  maximum  flow  network  be  called  the  primal.  To 
the  primal  add  an  artificial  arc  extending  from  the  source  to  the  sink 
and  having  a  capacity  of  zero.  The  resulting  network  will  be  referred 
to  as  the  modified  primal.  The  dual  network  is  defined  if  and  only  if 
the  primal  is  source-sink  planar,  a  source-sink  planar  network  being 
one  where  the  modified  primal  can  be  drawn  on  a  sphere  in  such  a  way 
that  no  two  arcs  intersect  except  at  a  node. 

When  defined,  the  dual  is  constructed  in  the  following  manner: 

1.  Draw  the  modified  primal  on  a  sphere  in  such  a  way  that  no 
two  arcs  intersect  except  at  a  node. 

2.  Place  a  node  in  each  mesh  of  the  modified  primal.  Let  the 
node  in  one  of  the  two  meshes  bounded  by  the  artificial  arc  be  the 
source  and  the  nodi  in  the  other  of  these  two  meshes  be  the  sink. 

3.  For  each  arc  except  the  artificial  one  construct  an  arc  of 
the  dual  t ha t  intersects  it  and  joins  the  nodes  in  the  meshes  on  either 
s ide  o  f  it. 

A.  Assign  each  arc  of  the  dual  a  length  equal  to  the  capacity  of 
the  primal  arc  it  intersects. 

Letting  a  route  through  the  dual  be  any  path  from  its  source  to  its 
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sink,  it  follows  that  there  is  a  one-to-one  correspondence  between 

[81 

the  proper  cuts  of  the  primal  and  the  routes  of  the  dual.  J  Specifically, 
if  A  is  any  route  through  the  dual,  then  the  arcs  of  the  primal  which 
intersect  A  form  a  proper  cut  and  conversely  if  B  is  a  proper  cut  of 
the  primal,  then  the  arcs  of  the  dual  which  intersect  B  form  a  route. 


-6- 


IV.  PROBLEM  FORMULATION 

As  stated  before,  the  arcs  of  the  maximum  flow  network  in  this 
paper  are  subject  to  anywhere  from  0  to  n  breakdowns.  The  decrease 
in  capacity  of  arc  (a,b)  due  to  i  breakdowns  is  a  random  variable 
with  unknown  distribution  but  with  known  mean  and  variance.  The  mean 

and  variance  of  this  quantity  will  be  denoted  as  u(a,b,i)  and 

2 

O  (a,b,i)  respectively.  Furthermore,  these  distributions  are  indepen¬ 
dent  for  the  different  arcs.  It  is  desired  to  find  the  smallest  possibl 
value  of  F  resulting  from,  at  most,  n  breakdowns  where  F  satisfies 

P{max  flow  >  F)  S0 

This  problem  can  also  be  formulated  in  terms  of  the  dual  network. 

Specifically,  each  arc  of  the  dual  is  subject  to  anywhere  from  0  to 

n  improvements.  The  decrease  in  length  of  arc  (a  ,b)  due  to  i 

improvements  is  a  random  variable  with  mean  u(a,b,i)  and  variance 
2 , 

O  (a,b,i)  .  It  is  desired  to  find  the  smallest  possible  value  of  F 
resulting  from,  at  most,  n  total  improvements  on  n  or  less  arcs 
where  F  satisfies: 

P{min  route  S  F}  ^  6 

It  is  the  latter  formulation  that  is  considered  in  this  paper. 

Since  examples  can  be  constructed  to  show  that  different  distributions 

2 

with  the  same  u(a,b,i)  and  a  (a,b,i)  can  have  different  F  and 
different  locations  for  the  improvements,  the  information  given  is  not 
sufficient  to  determine  a  solution.  Accordingly,  the  distribution 
considered  will  be  the  one  whose  smallest  value  of  F  is  maximum. 

Suppose  that  if  improvements  occur  on  certain  arcs,  the  length  of 
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a  particular  route  is  a  random  variable  L  with  mean  u  and  variance 

2  * 

n  .  It  follows  from  Tchebyshef f ' s  extended  lemma 

2 

p{l  -  u  ^  c}  <  -y2— 7 

c  +  a 


or 

2 

P{L  >  a  +  e}  <  —  ?  — 2 
e  +  a 

If  it  is  required  to  minimize  u  +  e  under  the  condition  that: 

tA‘. 

e  +  a 


it  follows  that 


and  one  obtains 


P(L  2  u  +  <p\J  y  1  }  <:  B 

Thu^  the  quantity  u  +  V  7:  -  1  plays  the  role  of  the  length 


of  the  above  mentioned  route  and  will  be  referred  to  as  its  effective 
length.  The  algorithm  to  be  presented  in  this  chapter  finds  those  arcs 
upon  which  the  n  improvements  must  occur  in  order  that  the  minimum 
effective  length  of  all  routes  through  the  dual  be  minimized. 

In  order  that  this  be  a  valid  criteria,  the  following  assumptions 
are  made : 


1.  u(a,b,i) 

2.  u(a,b,i) 


l(a,b)  and  is  strictly  increasing  in 

cr(a 


*Re f .  9,  pp.  111-126. 


1  >  0  for  i  i  1  . 
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3.  u(a,b,l)  >u(a,b,i)  implies 

u(a,b,i)  -  ©(a.b,!)^^  -  1  >  u(a 

These  assumptions  assure  that  the  solution  to  this  problem  involves 
exactly  n  breakdowns,  assure  that  the  effective  length  of  an  arc  is 
always  greater  than  zero  and  generate  arguments  in  favor  of  the 
efficiency  of  the  algorithm. 
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V.  THE  ALGORITHM 


Let  an  i-arc  path  to  node  a  be  any  path  from  the  source  to 
node  a  with  i  or  fewer  improvements  occurring  on  its  arcs.  Each 
node  is  assigned  n  +  1  sets  of  labels,  the  set  numbers  being 
designated  as  0,1 . n  .  Each  label  consists  of  four  components 

2  J 

and  is  denoted  (u,a  ,t,k)  .  where  a  designates  the  node,  i  the 

a  ,  i 

set  numb er,  and  j  the  rank  within  the  set.  The  individual  components 

are  denoted  as  ^  »  (cr^ )  ^  ,  t^  .  ,  and  k^  .  The  quantities 

d  1 1  3,1  a  ,  i  a|X 

1  2  1 

U  and  (ct  )  .  are  related  to  the  mean  and  variance  of  i-arc 

a  |  i  a  1 1 

1  2 

paths  to  node  a  and  at  termination  and  (o  )■*  are  the  mean 

)  ,n  t)  ,n 

and  variance  of  the  length  of  the  desired  path  where  j  ■  max  j  in  this 

set.  The  components  t^  .  and  k^  .  are  tracers  which  are  used  to 

a  ,  i  a  ,  i 

find  the  path  itself.  Again,  S  and  S  refer  to  the  source  and  sink 
respectively. 

The  algorithm  for  finding  the  mean  variance  of  the  length  of  the 
desired  route  is  as  follows: 


1.  Set  i  -  0  . 


2.  Set  ut  -  (<rV  -  0  .  ,  -  (<rV  . 

0,1  0,1  B  1 1  B  |  1 


for  a  i  S  and  s  ■  i 


3.  Consider  each  a  and  all  arcs  of  the  fora  (a,b)  .  For  each  j 
consider  the  quantity: 

L  "  Ub,i-s  +  1(a’b)  *  »*<afb,s) 


+  V(a2)b,i-s  +  a2(a'b’s)  "Vi 


1 
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and  delete  from  the  set  i  at  node  a  all  labels  j '  which 
satisfy: 


U 


U 


j' 

a  ,  i 


J* 
a  ,i 


-  u(a,b,s) 


L 


with  strict  inequality  for  at  least  one  of  these.  If  there  is  no  j' 
satis  fy  ing: 


i  €  i  c  +  l(a ’b>  "  u(a  *b«s) 

a , i  b  ,  l-s 


■  -j,.VTT 

,1  a,i  V  0 


£  L 


introduce  the  label  (fi  ,0  ,t ,k)  where 

a  ,  1 


“a,i  '  “b.t-s  +  1(a’b)  ‘  u(a'b's) 


Ccr2)a,i  =  ^b.i-s  +  ^ (a  'b  >s) 


t  .  -  b 


k  .  =  s 
a  ,i 

If  any  changes  result  in  the  labels  after  examining  an  arc,  re-rank  the 
labels  in  order  of  increasing  first  component.  This  results  in  a  change 
in  the  subscript  j  for  some  labels.  In  this  set,  j  now  takes  on  the 
values  l,2,...,r  where  r  is  the  number  of  labels  in  set  i  at  nodi  a. 

4.  If  s  =  0 ,  go  to  5  .  Otherwise,  decrease  s  bv  !  and  go  back 
to  3  . 

3.  Repeat  3  until  no  changes  in  the  labels  of  s*  t  i  at  node  a 


-11- 


result  . 

6.  If  i  n  ,  increase  i  by  1  and  go  back  to  2  .  If  i  3  n  , 

1  ,  2  1 

terminate,  as  ut  and  (a  )—  are  the  mean  and  variance  of  the  length 

^ i n  S  ,n 

of  the  desired  path  where  j  equals  the  number  of  members  of  set  n  at  S 
The  path  itself  may  be  found  by  the  following  procedure: 

1.  Set  m  3  1  ,  a  ^  =  S  ,  and  ij  =  n  . 

2 

2.  Let  (u,<7  ,  t  , k )—  '  1  be  the  label  of  highest  rank  among  set  n 

S  ,n 

at  S  . 

3.  Set  k  3  k^m 

m  a  ,  l 
m  m 

4.  Increase  m  by  1  .  Set  i  =  i  .  -  k  and  a  =  t^m"^ 

m  m-I  m-1  m  a  , , i 

m-1  m-1 

If  a  3  S  ,  go  to  6  .  Otherwise,  find  a  label  in  set  i  at  node  a 
m  mm 


which  satisfies 


*  - 


A  Jm 

3m*i 

m  m 


=  U 


Jm-1  .  -  l(a  .  ,a  )  +  U(a  ,a  ,k  ) 

a  ,  ,  l  ,  m- 1  m  m- 1  m  m- 1 

m-1  m-1 


(oVm  . 
VXm 


(a2) ^m-1  -  a2(a  ,a  ,k 

a  ,  ,  i  ,  m-1  m  nr 

m-1  m-1 


5 .  Go  back  to  3  . 

6.  Terminate.  S  3  a  . a,  3  S  with  arc  (a  ,a  )  having  k 

mi  J  J+i  j 

improvements  is  the  desired  path.  The  arcs  of  the  primal  upon  which 
breakdowns  must  occur,  are  those  which  intersect  this  path.  Specifically, 
k  breakdowns  must  occur  on  the  primal  arc  which  intersects  arc 

j 

(a^,a^+^)  in  this  path. 

*  2 

Note  that  u(a,b,0)  ■  c  (a,b,0)  ■  0, 


-12- 


VI.  JUSTIFICATION  OF  THE  ALGORITHM 


The  justification  will  consist  in  showing  that  the  algorithm  is 
finite  and  that  at  termination  the  label  of  highest  rank  in  set  n  at 
S  has  as  its  first  two  components  the  mean  and  variance  of  the  length 
of  the  n-arc  path  to  the  sink  of  minimum  effective  length.  Furthermore, 
the  steps  of  the  tracing  procedure  can  be  carried  out  and  finds  this 
pa  th . 


Lemma  1:  Let  S  ■  a„,a . .  with  k.  breakdowns  occurring  on  arc 

- —  0  1  m  j 

(a  ,a  )  be  an  i-arc  path  to  a  such  that  all  of  its  arcs  are  distinct. 
j-1  j  m 

Let  the  mean  and  variance  of  the  length  of  arc  (a^  ^,a^)  be  and 

2 

CTj  respectively.  Then  the  total  length  of  this  path  is  a  random  variable 


m 


m 


whose  mean  is  equal  to  Eu.  and  whose  variance  is  equal  to  E  cr  • 

j-1  J  j-1  J 

Proof :  The  amounts  by  which  the  lengths  of  the  different  arcs 

can  be  reduced  are  independent.  Hence,  the  lengths  of  the  arcs  themselves 
are  independent  and  their  variances  may  be  added  to  give  the  variance  of 
the  total  length  of  the  path.  Of  course,  the  means  may  be  added  regardless 
of  whether  or  not  the  distributions  are  independent. 


Corollary  2:  Let  S  *  a^,...  ,am  with  k;  breakdowns  occuriing 


m 


j 


on  arc  (a  ,a  )  be  a  non-cyclic  i-arc  path  to  node  a  ,  and  let  the 
J- 1  J  m 

2 

length  of  arc  (a^  paj)  have  mean  Uj  and  variance  .  Then  the 


m 


m 


length  of  this  path  has  mean  L  u  and  variance  Ecr 

j-1  J  j-1  j 


Prop  f :  Since  the  path  is  non-cyclic,  all  its  arcs  are  distinct  and 


the  corollary  follows  from  lemma  1. 
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2 

De  f  In  it  Ion :  The  pair  (u,a  )  is  said  to  dominate  the  pair 

- 2 

(li  ,n  )  if  and  only  if: 

i)  u  *  u 

if)  U  +  * 


with  strict  inequality  for  at  least  one  of  these.  If  path  1  is  an 

2 

i-arc  path  to  node  a  whose  length  has  mean  and  variance  , 


and  path  2  is  an  i-arc  path  to  node  a  whose  length  has  mean  u.^ 

2 

and  variance  •  then  path  1  is  said  to  dominate  path  2  if  and 


2  2 

only  if  the  pair  (u^  ,a^)  dominates  •  An  i-arc  path  to  node 


a  that  is  not  dominated  by  any  other  i-arc  path  to  node  a  is  said  to 

2  \ 

be  an  undominated  i-arc  path.  Furthermore,  a  label  (u ,ct  ,  t,k)  . 

a  ,  l 

2  1 

dominates  the  label  (u  ,o  ,t,k)  .  if  and  only  if  the  pair 

a ,  l 

6a\t  '  (°2>i,i>  d0”lnates  <5*1,1  ’  (<,2>],l)  • 

Note  that  the  dominance  property  is  transitive,  that  a  label  can 
be  introduced  into  a  set  only  if  no  other  label  in  that  set  dominates 
it;  and,  rinally,  when  a  label  is  being  considered  for  introduction 
into  a  set,  the  labels  that  are  dropped  from  that  set  are  precisely 
those  which  are  dominated  by  the  one  being  considered. 


Lemma  3 :  A  label  can  not  dominate  another  label  in  its  set. 
Prop  f :  Suppose  the  theorem  holds  at  one  stage  of  the  algorithm 
and  that  node  a  and  arc  (a,b)  with  s  improvements  are  being 

2  j 

considered  with  respect  to  the  label  (|i  ,c  ,t,k). 

D  *  1 "  8 
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Lt  t  : 

U  =  lij  +  l(a,b)  -  u(a  ,b  ,s) 

b  ,  l  -s 

a 2  =  (o2)^  .  +  a2(a  ,b  ,s) 

b  ,  l  -  s 

All  labels  among  set  i  at  node  a  whose  first  two  components  are 

- 2 

dominated  by  (U,0  )  are  dropped  from  this  set  and  the  condition 

2  1  - 2 

still  holds.  Then  the  label  (n  ,a  ,t,k)  =  (u  ,C  ,b,s)  is 

a  ,  i 

introduced  into  this  set  if  and  only  if  it  is  not  dominated  by  any 
other  label  in  this  set  and  the  property  is  still  preserved.  The 
only  other  way  a  change  in  labels  can  result  is  through  the  introduction 
of  (0,0,-,-)  or  C®, into  an  empty  set.  Initially,  the 
conditions  of  the  lemma  hold  with  all  sets  empty  and  the  lemma  follows 
from  induction. 

Lemma  4:  A  label  can  only  be  dropped  from  a  set  if  another  label 

which  dominates  it  is  introduced  into  that  set. 

2  j 

Proof:  Suppose  a  label  ,cr  ,t,k)  .  is  to  be  dropped  from 

a,i 

set  i  at  node  a  .  Then  there  is  a  quadruple  (p,,cr2,t,k)  which 
will  be  introduced  as  a  label  into  this  set  provided  no  other  label 
would  dominate  it  and  which  has  the  property  that  it  dominates 
2  1 

,a  ,t,k)  .  Furthermore,  no  other  label  in  this  set  can  dominate 

a  ,  i 

2  j 

it  for  if  it  did,  it  would  also  dominate  (u  ,o  ,t,k)  .  Hence, 

a  ,  i 

2  1  - 2  -  — 

(u  ,a  ,t,k)  .  =  (u ,o  ,t,k)  is  introduced  into  set  i  at  node  a 


proving  the  theorem. 


)  - 


Co  ro  1 1 .1  ry  r> :  If  a  label  is  dropped  from  a  p.irt  icn  l.i  r  sit  ,  there 
will  always  bo  a  label  in  that  sot  which  dominates  it. 

Prop  f :  Tit  i  s  follows  immediately  from  lemma  and  the  transitivitv 
of  the  dominance  property. 


Lemma  6:  Let  path  l,  S  =  an,...,a  =  a  with  k  improvi  mt  nt  s 

-  0  m  j 


on  i  .  c 

and  path  2  ,  S 

=  a„  ,  . .  .  ,a  =  a  with  k  . 
1)  m  ] 

imp rovenu  nt  s 

on  arc 

be  non-cyclic  i 

-arc  paths  to  node  a 

Li  t  path  1 

and  2  1 

be  (i+r)- 

arc  paths  to  node 

b  formed  by  adding  arc 

(a ,h )  with 

r  improvements  to  path  1  and  path  2  respectively.  If  path  1  dominates 
path  2  ,  then  path  1’  dominates  path  2'  . 

Proof :  Since  paths  1  and  2  are  non-cyclic.  it  follows  that  the 

arcs  of  each  path  1,  2,  1'  ,  and  2'  are  distinct.  Lett  in  }t  the  pairs 

(u(a)  ,a“(a))  ,  (u(a)  ,a2(a))  ,  (ji(b)  ,o2(b))  and  (jl(b)  ,rr"(b))  be  the 
means  and  variances  of  the  lengths  of  paths  1,  2,  1'  and  2'  ,  it 

follows  from  lemma  1  that: 

u(b)  =  u(a)  +  1  (a  ,b)  -  u(a,b,r) 

U(b)  =  u(a)  +  1  ( a  ,b)  -  u(a,b,r) 

2  2  2 

a  (b)  =  a  (a)  +a  (a,b,r) 

-2  -2  2 

a  (b)  =  a  (a)  +  a  (a  ,b  ,r) 

Therefore : 

u(a)  <  u(a) 

U(a)  +  1  (a  ,b)  -  u(a,b,r)  <  u(a)  +  l(a,b)  -  u(a,b,r) 


H(b)  <  n(b) 
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with  equality  only  if  u(a)  *  u(a)  and  the  first  condition  is  satisfied. 
For  the  second  condition  one  has: 

Case  1:  — ^  ____ 

u(a)  +  a(a)^jp-"T  =  u(a)  +  a(a)*^-g  -  1 

u(a)  <  u(a ) 
a(a)  >  a(a) 

2  —2  22  —22 
a  (a)  -  ct  (a)  *  [a  (a)  +o  (a,b,r)]  -  [a  (a)  +  cr  (a,b,r)] 

o^(a)  -  o^( a)  *  cr^(b)  -  a^(b) 
a(a)  +  cr(a)  <  a(b)  +  a(b) 
a(a)  -  a(a)  2*  a(b)  -  ar(b) 

U(a)  +  [a(a)  -  o(a)]^-g  -  1  =  u(a) 
u(a)  +  [a(b)  -  a(b)]*^A  -  1  S  u(a) 

U(b)  +  a(b T\j±  -  1  «;  u(b)  +  a(b)^  -  1 

Case  2: 

U(a)  +  a(a)^^  -  1  ^  ^i(a)  -t-  a(a)^|  -  1 

Define  a  such  that 
U(a)  + 

U(b)  + 

Also  o(a)  <  a  .  Thus: 

a(b)  <^/a2  +  a2(a  ,b  ,r) 


-17- 


li(b)  +  a(bf\/|  -  1  <  u(b)  +  a(b) 


VI 


Lenina  7:  For  each  label  (u,o^,t,k)^  ,  ,  the  quantities  , 

'  a  *  i  a » i 

2  \ 

and  (a  ) J  are  either  infinite  or  are  equal  to  the  mean  and  variance 
a  ,  l 

of  an  i-arc  path  to  node  a  . 

Proo f :  Suppose  at  one  atage  of  the  algorithm  that  for  each  label 
2  1  1  2  i 

(u,a  ,t,k)  where  u  ,  and  (ct  )  are  not  infinite,  there  is  an 

3,1  3,1  3,1 


i-arc  path  S  *  a  , . . .  ,a  =  a  with  k.  improvements  on  arc  (a.  ,  ,a  . ) 

0  m  j  J-l  J 

whose  arcs  are  distinct  and  whose  length  has  mean  and  variance 

a  ,  i 


2J 


m 


.  and  such  that  set  [i  -  £  k  ,1  at  node  a  contains  a  label 
a-‘  j=n+l  1 


whose  first  two  components  equals  or  dominates  the  mean  and  variance  of 

the  length  of  the  path  S  *  a a  with  k,  breakdowns  on  arc 

O  n  j 

2  "j 

(a^  i,aj)  •  Suppose  the  label  (u  ,a  ,  t,k)^  is  introduced  through 

the  examination  of  arc  (a,b)  with  r  improvements.  Then  there  is  a 
2  \ 

label  (u ,ct  ,t,k)  such  that: 

a  ,  i 

ub,l+r  *Ui,i  +  Ua'b)  '  “(a'b'r) 


<t,2)b,l+r  *  (<,2)i,I  +  a2(a’b'r)  ' 


Let  S  *  a 


0’’ ’  •  ’  m 


a  with  k.  improvements  on  arc  (a  .,a  )  be  a 


path  which  satisfies  the  above  conditions  for  (u  ,0  ,t  ,k) 


j-r  j 

j 

a  ,  i 


Cons ider 


the  path  S  -  aQ . a^ 


b  with  k  improvements  on  arc  (a  ,a  )  , 

J  J  J 


k  a  r  .  If  each  of  its  arcs  is  distinct,  it  satisfies  the  above 
m+1  * 
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2  \ 

conditions  for  (u  ,n  ,t,k)^  ,,  .  Suppose  its  arcs  are  not  all  distinct. 

b  ,  l+r 

Then  there  is  an  n  <  m  such  that  a  =  a  and  a  =  a  .  The 

n  m  n+1  m+1 

mean  and  variance  of  the  length  of  the  path  S  =  a^,...  >an+j  with 

m+1 

improvements  on  arc  (a  ,a.)  for  j  *"  n  and  V  k,  improvements  on 

j_1  J  j=n+l  J 

1  2  1 

arc  (a  ,a  )  dominates  the  pair  (Ul  . ,  >  (o  .  )  .  Also  the 

n  n+1  b,i+r  b ,  l+r 

path  S  =  a a  =  a  with  k.  improvements  on  arc  (i  ,a  ) 

On  j  J-l  J 

dominates  the  path  S  =  a  . a  with  k,  improvements  on  arc 

urn  j 

m 

(a.  ,a  )  •  Therefore,  if  £  k.  =  0  ,  there  is  a  labe  in  set  i  at 
J_1  j  j  =n+l  J 

node  a  which  dominates  (n,a^,t,k)^  .  contradicting  lemma  3  .  On  the 

a  ,  l 

m 

other  hand,  suppose  Z  k.  ^  0  .  Consider  any  label  in  set 

j=n+l  J 

m 

[i  -  Z  k .  J  at  node  a  whose  first  two  components  dominate  the  mean 
j=n+l  J 

and  variance  of  the  length  of  S  =  a  , .  .  .  ,a  =  a  with  k,  improvements 

On  j 

m+1 

onarc  (a  , a.).  This  label  and  arc  (a  ,b)  with  £  k.>r 

^  ^  j=n+l  ^ 

2  * 

improvements  will  be  examined  prior  to  the  introduction  of  (u>o  ,t,k)^ 

as  a  label.  After  this  examaination  there  must  be  a  label  among  set  (i+r) 

2 

at  node  a  =  b  which  equals  or  dominates  (u  ,a  ,t,k)^  .  •  But  this 

m+1  b  ,  l+r 


prevents  its  introduction  as  a  label.  Thus,  the  arcs  of  S  =  a  , . . . ,a  =  b 

0  m+1 


art  distinct.  The  introduction  of  (0,0,-,-)  into  any  set  at  the  source 
or  C00  ,  00  ,-,-)  into  any  set  preserves  the  above  properties.  Furthermore,  the 
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droppin^  of  a  label  preserves  these  properties  since  it  is  immediately 
followed  by  the  introduction  of  a  label  into  its  set  which  dominates 
it.  Initially  these  properties  hold  with  all  sets  of  labels  beinp, 
empty.  The  lemma  follows  frtxn  induction. 

Lemma  8 :  The  labeling  algorithm  is  finite. 

Prop f :  Since  the  first  two  components  of  a  label  are  either 
infinite  or  are  equal  to  the  mean  and  variance  of  a  path  whose  arcs 
are  distinct,  it  follows  that  this  pair  must  be  selected  from  a  finite- 

set.  In  addition,  the  choices  for  the  t^  and  .  are  finite. 

a  ,  1  a  ,  l 

Thus,  the  labels  themselves  are  selected  from  a  finite  set.  Furthermore, 
no  label  may  be  introduced  into  the  same  set  more  than  once,  for  if  it  is 
once  dropped  there  is  another  label  in  the  set  which  dominates  it  and 
prevents  its  re-entry.  Thus,  the  algorithm  is  finite. 

Lemma  9 :  At  termination  no  labels  are  infinite  at  any  node  to 
which  a  path  exists. 

Proof:  Consider  set  i  at  node  a  .  Let  S  =  a^ . a  =  a  be 

-  0  m 

s  t 

any  path  to  a  .  It  follows  that  on  or  before  the  (m+l)“™  examination 
of  the  nodes  for  set  i  in  step  5  a  finite  label  will  be  introduced 
into  the  set  i  at  node  a  . 

Corollary  10:  At  termination  .  and  (a^)^  correspond 

. . .  3  >  t  a  ,  l 

to  the  mean  and  variance  of  i-arc  paths  to  node  a  for  all  i  and  all 

j  provided  a  path  a  exists. 

Prop f :  This  follows  immediately  from  1  enrna  9. 

Lemma  1 1  :  Let  S  =  a^ . =  a  with  improvements  on  arc 
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(a^  i»aj^  an  undominated  i-arc  path  to  node  a  .  Then  this  path 
contains  no  cycles. 

P roo f :  Suppose  this  path  does  contain  a  cycle.  Let  (a^  ^,a^), 

J  ■  1 , . . . ,m  be  the  distinct  arcs  of  this  path.  Furthermore,  let  arc 

(a  ^  ,a  )  used  C  times  in  this  path  and  let  its  length  have 

J  J  j 

2 

mean  and  variance  a  .  Then  the  total  length  of  this  path  has 


m 


m 


.2  2 


mean  £  C  li  and  variance  £  C  ff  ,  C  ^  1  all  j  .  Deleteall 
J-l  J  J  J-l  J  j  j 

cycles  from  this  path.  The  resulting  path  is  an  i-arc  path  to  node  a 


m  _  m  _2  2  _ 

whose  length  has  mean  £  C  u  and  variance  £  C  cr  *  C  <  1  all  j 

j-l  J  J  j-l  J  J  J 


and  C  “  0  for  some  j  .  Thus,  the  path  S  -an,...,a  =a  is  not 


j 


m 


undominated. 


Lemma  12:  Let  S  =  a^,...,a  =  a  with  k,  improvements  on  arc 

—————  U  m  j 

(a  ,  ,a  )  be  an  undominated  i-arc  path  to  node  a  .  Then  S  -  a  , . . .  ,a 
j- 1  j  U  m- 1 

with  k,  improvements  on  arc  (a,  ,  ,a  )  is  an  undominated  (i-k  )-arc 

j  j-l  j  m 


path  to  node  a  . 

m-  i 

P roc  f :  Suppose  the  above  path  is  not  undominated.  Let 

S  ■  b_,...,b—  -  a  .  with  k  improvements  on  arc  (b  .  ,b  )  be  an 
U  m  m- 1  j  J - 1  J 

undominated  (i-k  )-arc  path  to  node  a  that  dominates  S  -  a_,...,a 

m  m- 1  U  m- 1 

with  k.  improvements  on  arc  (a.  .  ,a  )  .  Then  S  -  b  , . .  .  ,b  and 
j  J"1  j  0  m 

S  -  a„,...,a  ,  contain  no  cycles  and  it  follows  from  lemma  6  that 

U  m- 1 


S  -  bn,...,b_,a  with  k.  improvements  on  arc  (b.  ,,b.)  and  k 
u  mm  jr  j-l  J  m 
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improvements  on  arc  (b_ ,a  )  is  an  i-arc  path  to  node  a  which 

mm  m 

dominates  S  -  an,...,a  with  k.  improvements  on  arc  (a  ,,a.) 

0  m  j  j-1  j 

contradicting  the  hypothesis. 

I  2  1 

Lemma  13:  At  termination  the  set  of  pairs  .  ,(a  )J  is 

-  a*i  *»i' 

the  set  of  means  and  variances  of  all  undominated  i-arc  paths  to  node 

a  provided  a  path  to  a  exists. 

P roo f :  Suppose  that  after  the  algorithm  terminates  for  the  sets 

i  2  1 

r  the  lemma  holds  for  all  sets  of  pairs  (ii J  .  ,(o  ) J  ,  i  <  r  .  Let 

'  a  ,  i  a  ,r 

S  ■  a. . a  ■  a  with  k  improvements  on  arc  (a,  ,a  )  be  an 

0  m  j  j-1  j 

undominated  (r+l)-arc  path  to  node  a  and  let  (a  ,  ,a  )  be  the  last 

p-1  p 

arc  in  this  path  upon  which  improvements  occur.  From  repeated  application 
of  lemma  12  it  follows  that  S  ■  aQ,...,a^  ^  with  k^  improvements  on 

arc  (a.  .,8,)  is  an  undominated  (r+l-k  )-arc  path  to  node  a  . 
j-1  j  p  p-1 

Let  the  mean  and  variance  of  the  lengths  of  these  two  paths  be  the  pairs 

2  2 

(U  ,er  )  and  (u  ,  ,)  respectively.  Then  there  exists  a  pair 

mm  p-1  p-1 

(ui  , r+l-k  •('’2)i  ,  .r+l-k  >  ‘  Vl’ap-1>  '  th"'  ‘S  “  Ubel 

P-1  P  P-i  P  ^  P 

in  set  r+l-k  at  node  a  .  whose  first  two  components  are  u  and 

p  p-1  p-1 

2 

o  .  •)  Note  that  the  path  S  ■  a. . a  with  k±  improvements  on 

p-l  0  P  j 

arc  (a  ,1*.)  is  an  undominated  (r+l)-arc  path  to  node  a  .  Denote 
J" 1  J  P 

2 

the  mean  and  variance  of  the  length  of  this  path  by  u  and  O  .  There 

P  P 

j  2  j  2 

can  be  no  pair  (u  ,(a  )  . . )  that  dominates  (u  ,a  )  since  such 

a  *  rt*i  a  1  p  p 

p  p 
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a  pair  must  be  infinite  or  be  equal  to  the  mean  and  variance  of  an 
(r+l)-arc  path  to  node  a  .  Thus,  after  examining  set  (r+l-k  )  at 


node  a  ,  and  arc  (a  ,  ,a  )  with 

p-1  p-1  p 


2  j 

a  label  (u ,a  ,t,k)  f+1  such  t ha t 

3P’ 


k  improvements,  there  must  be 
P 


U 


1 

a  ,  r+1 
P 


2,  i 


U„  and  ("  'a  ,  r+1 


=  n 


Furthermore,  this  label  can  never  be  dropped.  Repeated  application  of 
this  argument  yields  the  conclusion  that  eventually  a  label 

(u  ,a  ,t,k)^  ,  will  be  introduced  where  u'  ,,  =  u  and  (n^)^  =  a" 

a,  r+1  a,  r+1  m  a,  r+1  m 

and  that  it  will  never  be  dropped.  Furthermore,  any  label  in  set  r+1  at 

node  a  whose  first  two  components  correspond  to  the  mean  and  variance'  of 

a  dominated  (r+1) -arc  path  to  node  a  will  be  dropped  bv  the  introduction 

of  a  label  whose  first  two  components  correspi  nd  to  the  mean  and  variance 

of  the  length  of  a  path  that  dominates  the  above  mentioned  one.  Thus  , 

when  the  algorithm  terminates  for  the  sets  r+1  ,  the  lemma  holds  for  all 

sets  of  pairs  (ll^  .  ,((7^)  ^  .)  ,  i  <  r+1  .  Let  S  *  b_ . b_  =  a  be  an 

r  N  a ,1  a  ,  v  0  m 

undominated  0-arc  path  to  node  a  .  An  argument  similar  to  that  above 


shows  that  after,  at  most,  m  +  1  checks  of  the  nodes  there  is  a  label 

(u  ,n  ,t  ,k)^  _  such  that  and  (a^)^  n  are  the  mean  and  variance 

a  ,  U  a  ,  u  a  ,  U 

of  the  length  of  this  path;  and,  finally,  when  the  algorithm  terminates 

j  2  j 

for  the  sets  0  ,  the  lemma  holds  for  all  sets  of  pairs  n’(c7  K  r)  ■ 

3  y  U  3  j  IT 

The  lemma  itself  follows  from  induction. 


Theorem  14:  When  the  algorithm  terminal,  s,  the  label  of  highest 
rank  among  set  n  at  S  has  as  its  first  two  components  the  mean  and 
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variance  <  f  the-  n-arc  path  of  minimum  effective  length  to  S  . 

Prop  f :  The  n-arc  path  of  minimum  effective  length  is  undominated  and 

therefoic  there  must  exist  a  label  (a  ,a"  t  ,k)i-  at  termination  such 

S  ,n 

that  uj-  and  (rj^)  —  correspond  to  t he*  mean  and  variance  of  the 
S  ,  n  S  ,  n 

length  of  this  path.  Furthermore,  there  can  be  no  other  label 


9  * 

(li  ,a  ,t  ,k)^  of  higher  rank,  for  if  so  one  has: 
S  ,  n 


4 

5  ,  n 


4,n 


contradicting  lemma  3. 

Of  course  the  n-arc  path  of  minimum  effective  length  to  S  is 
the  route  ot  smallest  effective  length  possible  from  source  to  sink 
due  to,  at  most,  n  breakdowns.  Thus,  the  algorithm  does  find  the 
mean  and  variance  of  the  desired  path.  It  now  must  be  shown  that  the 
tracing  procedure  actually  finds  this  path  and  the  arcs  upon  which  the 
n  Improvements  must  occur. 


2  i 

Lemma  15:  Let  (u,a  ,t  ,k)  .  ,  a  i»  S  be  a  label  at  termination. 


2  j 

Then  there  is  another  label  (u  ,n  ,t,k),  —  at  termination  such  that: 

b ,  i 


b  *  t 


J 


a  ,  l 


i  *  i  -  k 


a  ,  i 


J 


\  '  4  7  +  l<a,b)  -  i*<a,b,k-j  ) 
a  ,  i  b  ,  i  a  ,  l 


■  Co'>bj  +  c’/(a>biki.i> 


2J 
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2  J 

Proof:  Consider  the  label  from  which  (u ,a  ,t,k)_  J  was  obtained. 

'  a  » » 

2  1 

This  label  must  satisfy  the  above  conditions,  i.e.,  must  be  (u ,a  ,t,k)^  — 

Furthermore,  this  labM  can  not  be  dropped,  for  if  it  <s  another  label 
2  1* 

(li,a  ,t,k)^  T  introduced  which  dominates  it.  Examining  this  label 
b  ,  i 

and  arc  (a,b)  with  k^  improvements  produces  a  label  which  dominates 

a  ,  l 


and  hence  drops  (u,a^,t,k)^  . 

a  ,  i 


Lemma  1 6 :  The  tracing  procedure  is  finite. 

Prop f :  The  sequence  of  labels  found  is  strictly  decreasing 

lexicographically  in  ,  i)  and  hence  none  can  be  repeated.  Thus, 

a  ,  i 

the  process  is  finite. 


Lenina  17:  Let  the  sequence  of  labels  found  by  the  tracing 
procedure  be: 


(U,cr  ,t  ,k)^r  ,  r  «  1 ,. . .  ,m 

a  |  1 

r  r 


Then: 


m-  1 


uli  -  u-h  -  \ ' 

S  »n  a^  ,i^ 

r-1 


[1(®r»ar+i^ 


u(V*r+l>kr)] 


m- 1 


}.i,  '  2^2<ar’ar+l-kiJ 


where  k  ■  k  r  , 
r  a  ,i 
r  r 


Proof:  It  follows  from  the  rules  of  the  tracing  procedure  that: 
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<‘>  “ar,l  '  “ir+!,l  x  ‘  1(V\+1>  •  “‘Vr+l’V-  r  '  ' . m' ' 

r  r  r+1  r+1 


m  ii  Jrn  ■  0 

S  ,  i  a  ,i 
m  ram 


r2\  -^r 


A  ^r+1 


(ii)  (a')^r  '  (a')pr+1  <  -  a  (a  ,a  k  ) ,  r  -  1,. . .  ,m-l 

ar’4r  ‘ r+1'1 r+1  r  r+1  r 


(CT2,^  -  (a2)J™  -0. 

ra  mm 


Sunning  the  equations  (i)  gives  the  first  relationship  and  summing  the 
equations  (ii)  gives  the  second. 

Theorem  18:  Let  the  sequence  of  labels  found  by  the  tracing 

procedure  be  as  in  lemma  17.  Then  S  *  a  , .  .  .  ,a  -  S  with 

ra  1 

k  ■  k^r  ,  improvements  on  arc  (a  ,a  , . )  is  an  n-arc  path  from 
r  a  ,i  r  r  r+1 

r  r 

source  to  sink  whose  length  has  mean  and  variance  (a 

S  ,n  S  ,n 


Proof:  The  quantities  and  equal  the  sums  of  the 

means  and  variances,  respec  t  ive  ly,  o  f  the  arcs  in  this  path.  All  the  arcs 
in  this  path  are  distinct,  for  if  not  there  would  be  a  cycle;  and, 
deleting  all  cycles  would  yield  an  n-arc  path  to  S  such  that  its 

length  has  mean  strictly  less  than  lir*  and  variance  not  exceeding 

S  ,  n 

2  i 

(a  contradicting  theorem  14.  Thus,  the  theorem  follows  from 

b  ,n 

leirma  1 . 

This  completes  the  justification  of  the  procedure  for  finding  the 
desired  route  itself.  The  primal  arcs  upon  which  breakdowns  must  occur 
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are  those  which  intersect  this  route.  In  particular,  k,  breakdowns 


must  occur  upon  that  arc  of  the  primal  which  intersects  arc 


‘Wi’ 
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VII.  NORMAL  CASL 


Suppose  it  is  known  that  the  decrease  in  the  capacities  of  the 

arcs  due  to  breakdowns  are  independently  and  normally  (or  rather 

approximately  normally)  distributed  with  known  mean  and  variance. 

It  follows  that  for  any  breakdown  pattern  all  routes  through  the 

dual  are  normally  distributed. 

Let  the  length  of  a  particular  path  be  a  random  variable  L 

2 

whose  distribution  is  normal  with  mean  u  and  variance  a 
Then : 


>  ♦(!-»)  j 


!  . 


where  $(1-8)  satisfies 


P  j Y  £  $(1-8)  J  -  1  -  8 

where  Y  is  standard  normal.  Thus: 


P  j  L  >  u  +  CT$(1  -  0)  J  =  8 

Hence,  u  +  ct$(l  -  8)  plays  the  role  of  the  effective  length  of  this 
path.  Suppose  it  is  desired  to  find  the  n-arc  path  to  the  sink  of 
minimum  effective  length.  Then  the  algorithm  just  presented  may  be 

used  to  solve  this  problem  by  replacing 


1  with  $(1  -  8) 
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VIII.  PRACTICALITY  OF  THE  ALGORITHM 

Since  this  algoritlim  is  in  many  ways  similar  to  the  minimum 

•k 

rcute  algorithm  which  is  highly  efficient,  it  would  seem  t ha t  this 
algorithm  would  be  practical,  provided  the  number  of  labels  in  each 
set  remained  small.  While  examples  have  not  been  solved  to  establish 
whether  or  not  the  number  of  labels  in  each  set  is  likely  to  remain 
small,  there  are  indications  in  favor  of  this  happening. 

The  first  of  these  indications  is  the  assumption  that: 

u(a  ,b  ,  i )  >  u(a  ,b  ,  i) 


implies 


U(a  ,b  ,  i )  -  o( a 


•b-i>Y!  ■ 


1  >  u(a  ,b,i)  -  o(  a  ,b  ,  i  ) 


Thus,  the  difference  between  the  zero-breakdown  length  of  a  path  and 
its  effective  length  tends  to  increase  as  the  mean  decrease  in  length 
due  to  breakdowns  increases. 

The  second  of  these  indications  is  even  more  convincing.  Let 
A  and  B  be  two  paths  to  node  a  such  that  A  has  smaller  mean 
length  than  B  but  does  not  dominate  it.  Then  if  arcs  with  breakdowns 
are  added  to  the  two  paths,  the  new  path  formed  from  A  tends  to 
dominate  that  formed  from  B  .  This  is  summarized  in  theorem  1. 


Theorem  1:  Consider  the  following  quantities: 


U 


1 


U 


2 


+  U 
U2  +  14 


★ 

Ref.  7,  pp. 


130-134. 
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2  2  —  —  —  — 

Then  if  ,  U j_  "  ^2  *  ul  “  ^2  and  L1  "  L2  ^  L1  ’  L2  with 

2 

equality  if  and  only  if  a  ■  0  . 

Proof : 


-  u 2  “  *  ^2  ^°^ows  trivially 

(L1  -  L2)  -  (l{  -  l2) 


2 

Ncte  that  equality  holds  in  the  above  if  a  ■  0  . 
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The  significance  in  'his  is  that  it  tends  to  keep  the  number  of 
labels  in  a  set  from  increasing  at  a  rapid  rate  as  one  gets  further 
from  the  source. 
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LIST  OF  SYMBOLS 


(a,b) 

Arc  joining  a  and  b 

1 (a  ,h) 

Length  of  arc  joining  a  and  b 

u(a  , b  ,  i ) 

Mean  of  the  decrease  in  capacity  of 
breakdowns 

(  a  ,b )  d ue  to  i 

2 

a  (a ,b  ,  i ) 

Variance  of  the  decrease  in  capacity 
b  reakdowns 

o  1  ( a  ,b)  due  to 

n  • 

(u  ,CT  ’  ,  t  ,k)  1 

a  ,  l 

A  label  for  the'  stochastic  case 

U  ’  . 

2 

First  component  of  (u  ,rr  ,t,k)' 

a  ,  i 

a  ,  i 

(a2)'  . 

2  j 

Second  component  of  (u,CT  ,t,k) 

a  ,  i 

a  ,  l 

tj  . 
a  ,i 

2 

Third  ce>mponent  of  (u  ,n  ,t,k) 

a  ,  i 

kj  . 
a  ,  l 

2  j 

Fourth  component  of  (u.ct  ,t,k) 

a  ,  i 

[A.B] 

Cut  set  [A,B] 

V[A,B] 

Value  o  f  cut  set  [A,B] 
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